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Physical contribution of the high-energy part of hadronic spectrum is incorporated 
to describe the short-distance part of the correlator of heavy-quark currents ob- 
tained by quenched Monte-Carlo on a 8 3 * 16- and 16 3 * 32-lattices for /3 = {6, 6.3}. 
The lattice artifacts in the short-distance behavior of that correlator are isolated. 
The physical short-distance part of the correlator is fitted by the relevant expres- 
sions of perturbative QCD, which allows one to obtain the renormalized charmed 
quark mass with rather high accuracy m^ s (m c ) — 1.22(5) GeV. 



1 Correlator of the heavy-quark currents 

The correlator of the vector currents of charmed quarks: 

n^V = ijdxe* x <0\T{j ll (x)j v (0)}\0> (1) 

where = cj^c can be used to calculate the parameters of perturbative 
QCD, such as the renormalized charmed-quark mass m c and the strong cou- 
pling constant a s (4m 2 ) 1 . On one hand this correlator satisfies the standard 
dispersion relation: 

ds ~( — I — 2T + dl 9v + dzilnlv -9nvQ 2 ) 
s ys + q ) 

(2) 

with the singular subtraction constants d\$,. The spectral density p(s) is pro- 
portional to the observable inclusive cross-section of charm-anticharm produc- 
tion in the e + e~-annihilation. It is well described by the following simple 
ansatz 1 : 

Pphen(s) = s (^fm 2 , e J(s - m 2 res ) + ^ d ( s - s o)^J ( 3 ) 

with the low- lying resonance of mass m res = mj/^, residue / proportional to 
the electromagnetic width of the J/^-meson, and a smooth continuum spec- 
trum with some effective threshold s > m 2 es . The experimental values of the 
spectrum parameters are: 

4tt 2 / ~ 0.6, m res ~ 3.1 GeV, s ~ (4 GeV) 2 (4) 
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On the other hand this correlator is computable in perturbative QCD in the 
vicinity of vanishing momentum q 2 = as that point is far away from the 
nearest threshold Am 2 due to charmed quark-antiquark pair. The parameter 
of the perturbative loop expansion here is a s (4m 2 ) < 1. Following 1 we 
consider the ratios r n = M n +i/M n of moments of the correlator (1). The 
moments are defined as 



The applicability of perturbative QCD near q 2 = implies the following ex- 
pansion for the ratios r n : 



where {a n , b n ,c n } are known numbers 1 . The term ~ a n comes from one 
loop of free charmed quarks. The term ~ b n comes from the two- loop dia- 
grams corresponding to one-gluon exchange. The gluon condensate = 
< 0|(a s /7r) G^„G^„|0 > is the vacuum average of the first nontrivial operator 
in the Wilson operator-product expansion for the correlator of heavy-quark 
currents. The coefficient c n originating from the Wilson coefficient function 
starts with one loop of the heavy-quark propagators. 

The lower ratios r2,3,4, originating from short distances, are well repro- 
duced in perturbative chromodynamics. The typical virtuality corresponding 
to the nth moment is ~ Am 2 /n. The coefficient b n grows with n. At large 
n > 8 the perturbation-theory based expansion (6) is irrelevant. The gluon 
condensate shows up in the intermediate ratios r^^^ 1 sensitive to larger dis- 
tances and, hence, nonperturbative fluctuations. The lower moments receive 
significant contribution from the high-energy (continuum) part of hadronic 
spectrum, whereas the higher moments correspond to long distances and are 
saturated by the lowest resonance J/tp. The moment Aio is quadratically 
divergent and renormalizes the photon mass, whereas the moment Aii is loga- 
rithmically divergent and renormalizes the photon wave function. All the other 
moments are physical quantities: their ultraviolet singularities are absorbed by 
the renormalization of the QCD Lagrangian parameters, the charmed-quark 
mass and the strong coupling constant. 

One can obtain a good estimate 1 of the renormalized QCD Lagrangian 
parameters by comparing the perturbative expressions for the lower ratios (6) 
to the same quantities obtained phenomenologically on the basis of dispersion 
relation (2) . The principal uncertainty of this estimate is that the continuum 
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contribution, significant in the lower ratios, is not known well enough from 
experiment. It is also only an approximation to identify the spectral density of 
the correlator of charmed-quark currents with the inclusive cross-section of the 
charm-anticharm production in e + e~-annihilation because of the contribution 
of light-quark currents to that inclusive cross-section. To avoid this problem 
and obtain an accurate way of calculating the renormalized charmed-quark 
mass and the strong coupling constant normalized on the charmed threshold it 
was suggested in 2 , 3 to calculate the correlator (1) by Monte-Carlo in Lattice 
QCD. Then the only experimental information needed is the low-lying reso- 
nance mass (mj/^). It is used to fix the scale. This idea can only work if one 
can isolate physical contribution of the high-energy part of hadronic spectrum 
in Monte-Carlo data on the background of short-distance lattice artifacts. We 
demonstrate below that this is indeed the case. 
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Figure 1: Ratios of the neighboring moments of the single resonance approximation to the 
correlator, on three lattices of different sizes: 8 4 , 16 4 , 32 4 . The resonance mass (m res = 0.5) 
is given in the units of the inverse lattice spacing. Solid lines are to guide the eye. Dashed 
line shows the moments of the subtracted correlator. Horizontal dotted line is a prediction 

of the continuum theory. 



2 Lattice artifacts of the moments 

We plan to analyse Monte-Carlo data obtained on relatively small and coarse 
lattices. Hence the lattice artifacts is the number-one problem to deal with. 
As usual we distinguish finite- volume from the cut-off effects. 
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Rather than doing Fourier transform of the correlator (1) and evaluating 
derivatives of (5) on a discrete lattice we calculate the moments directly in the 
X-space: 

We consider the correlator (1) in the inflnitly-narrow-resonance approximation 
to be the propagator of a free field of mass m res : 

n res (a) = - —r^-nn — - (8) 
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The moments of a single resonance (8) computed with the help of eqn.(7) on 
lattices of different size are shown on Fig. 1. The straight horizontal line is 
expected in the continuum theory r^ es = \jm 1 res . One can see strong finite- 
volume effects on small lattices 8 4 , 16 4 which persist even on the big (in a 
computational sense) 32 4 -lattice in the high moments. The higher moments 
originate in the long-distance part of the correlator. One can see the horizontal 
plateau of the continuum theory show up in the lower moments on the big 32 4 - 
lattice. On bigger lattices that plateau will expand to higher moments, but 
it will always remain somewhat above the prediction of the continuum theory 
due to cut-off 0(a 2 )-effects. Those cut-off effects in the moments of a single 
resonance are negligible for a ■ m res < 1 3 . 

The lower moments of the correlator (1) are expected to be well reproduced 
in the loop expansion of perturbative QCD, where the dominant term is the 
one loop of free quarks with the renormalized mass rh c . We need to know 
therefore the lattice artifacts of one loop of free Wilson fermions (as the Monte- 
Carlo data analysed was obtained for Wilson fermions). The moments of the 
correlator (1) in the one-loop (free quarks) approximation are shown in Fig. 
(2) for lattices of various sizes. Again one finds finite- volume effects which are 
particularly strong in high moments. Even on the big 32 4 -lattice the moments 
of one loop of Wilson quarks deviate from the continuum theory prediction 
(diamonds) due to the cut-off 0(a)-effects. One can reconcile nicely the lattice 
moments (for L > 32) with the expectations of the continuum theory by using 
only a simple relation between the quark mass of the continuum theory m c 
and the mass parameter of the Wilson propagator m c = 1/2/t — 4 3 : 

m 2 c = m 2 c I (I + m c ) (9) 

All the masses in (9) are in the units of a. The mass fh c is defined from the 
behavior of the quark propagator near vanishing momentum. It has nothing 
to do with the pole mass. 
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Figure 2: Ratios of the neighboring moments of the subtracted (dashed lines) and unsub- 
tracted (solid lines) correlator of pseudoscalar currents, saturated by one loop of free Wilson 
fermions with k = 0.1 (am c = 1) and periodic boundary conditions, computed numerically 
on lattices of indicated sizes. Lines connect the data to guide the eye. The diamonds is a 
prediction of the continuum theory with the quark mass m c = m c . 



Note that the relation (9) describing 0(a)-fermionic cut-off effects is valid 
in the infinite volume. Corrections to this formula due to finite volume were 
studied in 3 . 

Strong finite- volume effects imply that one has to fit Monte-Carlo data not 
with the continuum-theory expressions for the moments but with the moments 
of a single resonance or free quarks, calculated in the same box where the 
Monte-Carlo simulation was done, which is what we do. 



3 Dispersion relations on the Lattice 

Our aim is to fit the short-distance part of the correlator of heavy-quark cur- 
rents with the expressions of perturbative QCD. The short-distance part of a 
two-point correlator is determined to a great extent by the high-energy part 
of the spectrum. The problem therefore is to distinguish between the physical 
contribution of the high-energy part of spectrum and the pure lattice arti- 
facts which will eventually dominate as one moves into the domain of very 
short distances. We claim to see a window of intermediate distances where the 
contribution of the high-energy part of sectrum is significant while the lattice 
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artifacts do not dominate. Incorporating the high-energy part of the spectrum 
into the fit of Monte-Carlo data for the two-point correlators amounts to gen- 
eralization of the dispersion relation to the lattice theory. The problem, in 
other words, is to construct the correlator of the lattice theory corresponding 
to a given spectral density. 

We explore the following extention of the dispersion relation for the two- 
point correlator II(g 2 ) to the lattice theory: 

. 4 

U(q 2 ) = di + d 2 — ^sin 2 (^a/2) (10) 

+ f4y> 2 (^)V/^ — (ii) 

V, ^2 'J J ** (* + ^E>m 2 (W2)) V 
The factor ~ sin 2 is the inverse of the discrete Laplacian: 

(12) 

Therefore the coefficient di in (10) is precisely the renormalization of the pho- 
ton wave function in the case of vector currents of charmed quarks. The 
renormalized correlator would be given by eqn.(10), (11) with d\ = d% = 0, 
considered as a function of the renormalized Lagrangian parameters (we also 
call it "subtracted correlator"). The subtraction constants d\$, are ultraviolet 
singular. In contrast to the continuum theory the term against cfe in (10) is 
not a polynomial in q 2 , so it contributes to all the moments, not only to M\. 
We find that its contribution to M n decreases very fast with n. If the quark or 
resonance mass in units of a is less than 1, the ratio r% is affected significantly 
by the term ~ di (in contrast to the continuum theory, where this ratio is 
physical). The ratio is affected very little, whereas the contamination of all 
the other ratios r„>3 is negligible. One can see the difference between the sub- 
tracted and unsubtracted correlators on Figs. (1) and (2). The short-distance 
lattice artifacts are confined to the ratios 7*1,2 with rather small effect on as 
soon as the relevant masses in units of a are less than 1. We therefore choose 
to fit the ratios rs and with the expressions of perturbative QCD. 

4 Monte-Carlo data 

The ratios of the moments of the correlator of vector currents of the continuum 
theory 1 are shown on Fig. 3. Higher moments of the correlator, obtained on 
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Figure 3: Ratios of the neighboring moments in the continuum theory for the correlator of 
interpolating currents of the J/ip-meson. Burst symbols correspond to the phenomenological 
evaluation of the correlator via dispersion relation with the smooth hadronic continuum 
spectrum modelled as in [2]. Squares show the contribution of a single J/i/)-resonance. 
Diamonds are moments of one loop of free quarks of mass m c = 1.26 GeV. 



the basis of the dispersion relation and experimental information about the 
spectral density, approach the moments of a single resonance (J/^-meson), 
which always lie on a straight horisontal line. The higher moments therefore 
come from long distances. The lower moments (7*2,3 ) are seen to be fitted well 
by the moments of one loop of free charmed quarks, they are short-distance 
quantities. The two-loop correction due to one-gluon exchange will improve 
agreement between perturbative predictions and experimental predictions of 
the lower moments of the correlator. Incorporation of the gluon condensate 
will improve fit in somewhat higher moments 7*5,6,7 1 - 

Compare Fig. 3 with Fig. 4, which shows the ratios of the successive 
moments of the same correlator computed by Monte-Carlo on 16 3 ★ 32-lattice 
for f3 = 6, k = 0.1100 with the tadpole-improved Wilson Clover action 4 . 
One can see strong finite- volume effects in Monte-Carlo data in the domain of 
higher moments. The moments of a single-resonance (solid line) exhibit the 
same finite-volume effects and fit Monte-Carlo very well in a wide range of 
higher moments. This allows us to fix the scale of the lattice theory. 

The lower moments {r n ,n < 6} are seen to deviate strongly from the 
single-resonance curve. All Monte-Carlo moments are fitted remarkably well 



7 



0.6 



S-H 

0.4 



0.2 

5 10 

n 

Figure 4: Monte-Carlo data on the 16 3 x 32, /3 = 6 lattice for vector currents, k = 0.1100. 
The solid line is a single resonance approximation. The dashed line has been obtained with 
the phenomenological charmonium spectrum incorporated via the dispersion relation. 

by the dashed line. This line shows the moments of the lattice correlator 
corresponding to the spectral density which incorporates a smooth continuum 
spectrum in addition to the low-lying resonance. That lattice correlator is 
obtained using the disperion relation (11). The parameters of the spectral 
density are chosen to have experimental values as in (4). Although the lattice 
artifacts in Monte-Carlo data are strong, they are well described by the natural 
modification of the denominator in (11), which describes propagation of a state 
of a given invariant mass s on a finite lattice. We find evidence that it is a 
good approximation to ignore possible lattice artifacts in the spectral density, 
as it is done in (11). Incorporation of the high-energy (continuum) part of 
hadronic spectrum changes the line of ratios r n dramatically in small n. As 
expected the high-energy part of spectrum is important to descibe the short- 
distance behaviour of the correlator. It is less trivial, that one can identify the 
physical contibution of the continuum spectrum and distinguish it from the 
short-distance lattice artifacts, which are confined to the lowest ratios ri^, on 
the lattices of reasonable size. 

5 How to fit Monte-Carlo data at short times 

The successfull incorporation of the continuum spectrum helps to fit Monte- 
Carlo data on shorter distances and hence fix the scale on small lattices. Fig. 
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Figure 5: Resonance mass from the conventional fit of the zero-momentum component of 
the correlator of vector currents to cosh[(t — L/2)m re s] on the smaller 8 3 x 16 lattice. 

5 shows the conventional way of extracting the lowest-resonance mass from 
the two-point correlator of local currents on the small 8 3 * 16-lattice: Monte- 
Carlo data are fitted by the expression for the correlator in the signle-resonance 
approximation. Horizontal plateau is expected at large times which height is 
given by the resonance mass. That plateau is not seen on Fig. 5 because the 
8 3 *16 lattice is too small. The same type of picture holds (see Fig. 6) in terms 
of the moments: the resonance mass exptracted from the single-resonance fit 
to different moments should exhibit a horizontal plateau. On the 8 3 *16 lattice 
that plateau only starts to form at higher moments. At lower moments the 
resonance mass deviates from the horizontal line. This deviation is due to 
significant contribution of the high-energy part of hadronic spectrum at short 
distances. 

We now include the continuum spectrum as in eqn. (4) via the dispersion 
relation (11) keeping the ratio of the effective continuum threshold to the 
resonance mass fixed to have the experimental value s jm 1 res « 1.7 and use 
the resonance mass as the only fit parameter. The result of such a fit is shown 
on Fig. 6. One can see that the incorporation of the continuum spectrum via 
the dispersion relation (11) allows one to fix the resonance mass unambiguously 
from the data on relatively short distances. 

Proceed now to a bigger 16 3 * 32 lattice, where one can determine the 
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Figure 6: Effective resonance mass from fitting the moment ratios of Monte-Carlo data on 
the 8 3 x 16, j3 = 6 lattice for vector currents, k = 0.1060. The upper line is a single-resonance 

approximation. 

resonance mass by the conventional fit anumbiguously since the time extent is 
large enough for the high-energy excitations to decay. The ground state dom- 
inates and constitutes horizontal plateau at t > 10. The resulting resonance 
masses are shown in Tables 2, 3. One can compare them to the resonance 
masses extracted from the small lattice by means of the dispersion relation 
(incorporating the continuum spectrum), shown in Table 1. The agreement is 
very good. The corresponding spectrum of charmonium bound states is shown 
on Fig. 8. 

6 The renormalized charmed-quark mass 

We fit the lower moment ratios rz t i of Monte-Carlo data with the free-quarks 
approximation to the correlator (1): 



The Wilson parameter k corresponding to those free quarks determines the 
the renormalized quark mass in accordance with (9) if Monte-Carlo is done on 
sufficiently large lattice, such as 16 3 *32. On smaller lattices the relation (9) 
needs to be corrected due to finite-volume effects. This is done by means of 
the dispersion relation 3 . 
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Figure 7: Resonance mass from the conventional fit of the zero-momentum component of 
the correlator of vector currents to cosh[(t — L/2)m re s] on the bigger 16 3 x 32 lattice. 



In order to fix the subtraction scheme one should be sensitive to the a s - 
correction in the ratios r n . The study of the lattice artifacts on the two- loop 
level is not done yet. However, since the term b n a s of eqn. (6) is small in the 
lower ratios r^,A, its value may be taken from the continuum theory: 

- a " t 

£ n = 1 + a s b n /a n (14) 

We take the value of the strong coupling constant a s (m c ) « 0.3 and the co- 
efficients b n /a n - from the two- loop calculations of the continuum theory 1 for 
the vector, pseudoscalar and scalar channels: 

$ sc = 1.02 , $ ec = 0.95 , £ 3 sca = 0.97 

£f c = 0.96 , Q ec = 0.92 , £f ° = 0.91 (15) 

The values of £ n in (15) correspond to our choice of the coefficients b n in the 
M S'-subtraction scheme with the normalization point /x = m c 5 . After the fit 
of Monte-Carlo data with free fermions (one-loop approximation) the mass of 
those fermions must corrected by the coefficients £ in accordance with (15). 

The results for the renormalized charmed-quark mass are shown in Tables 
1, 2 , 3. The mass is extracted from the ratios rs and r± independently in 
three different channels: pseudoscalar, vector and scalar; on two lattices of 
size 8 3 * 16 and 16 3 * 32, for a few values of the hopping parameter k in the 
vicinity of n c harm of which k = 0.1060 (for (3 = 6) and k = 0.1150 (for (3 = 6.3) 
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Figure 8: Charmonium spectrum obtained on the 16 3 X 32 lattice at /3 = 6 (squares) and 
j3 = 6.3 (diamonds) and on the small lattice 8 3 x 16 at /3 = 6 (crosses) with the clover-and- 
tapole-improved action. The pseudoscalar mass m rlc — 2.979GeV is used for normalization. 
The wide horizontal lines are experimental data. 



are shown. The values of k correspond to the tadpole-improved Wilson Clover 
action 4 . The data at /? = 6.3 checks the scaling properties of the renormalized 
heavy-quark mass. By definition this physical mass should stay fixed and finite 
as a — > ({3 — > oo). 

One can see a rather stable value of the renormalized charmed-quark mass 
mf s (ffl c ) = 1.22(5)GeV, which is well in agreement with the estimates of the 
continuum theory 6 : m^ s (m c ) ss 1.23GeV. The previously reported lattice 
result is 7 : mf s (m c ) = 1.5(3)GeV. 

7 Conclusions 

Study of the short-distance behavior of the correlator of heavy-quark currents 
on the lattice appeares as a very promissing way to accurately calculate the 
parameters of perturbative QCD. 
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Table 1: The charm-quark mass obtained from fitting the ratios r^,, ta of the subtracted 
correlator on the small lattice 8 3 x 16 for f) = 6j k = 0.1060. The masses m re s and tn c are 
in units of the lattice spacing. The masses m^ ls [a J /^\ and m ( ^ f " s [a^ e ] are in GeV . They 
are obtained assuming that the lattice spacing is fixed from the vector channel and from the 
pseudoscalar channel respectively (a -1 i=» 1.9GeV). The indicated errors are statistical. 



fit tO T3 


7Tl res 


rh c 


Tile / "Wires 






Vc 
Xo 


1.57(1) 
1.61(2) 
1.87(10) 


.596(2) 
.651(4) 
.609(9) 


.381(2) 
.405(3) 
.329(10) 


1.16(1) 
1.22(1) 
1.15(3) 


1.15(1) 
1.21(1) 
1.14(3) 


fit tO T4 


7Tl res 


rh c 


TYt ^ 1 fy^if g § 


m^ s [aj^} 


m^ s [a nc ] 


Vc 

J/V 

Xo 


1.57(1) 
1.61(2) 
1.87(10) 


.623(3) 
.666(4) 
.649(10) 


.398(2) 
.415(3) 
.350(10) 


1.17(1) 
1.23(1) 
1.19(3) 


1.16(1) 
1.22(1) 
1.18(3) 


Table 2: Same as in Table 1 for the big lattice 16 3 * 32, /3 = 6, k = 0.1060. 


fit to rs 


7n res 


m c 


rh c /m res 




mf[«,J 


Vc 

j/4> 

Xo 


1.562(5) 
1.600(6) 
1.86(4) 


.618(2) 
.660(1) 
.645(2) 


.396(2) 
.413(2) 
.347(2) 


1.208(6) 
1.246(6) 
1.230(8) 


1.190(6) 
1.228(6) 
1.212(8) 


fit tO T4 


77l res 


m c 


th c J Tn res 




mf s [a Vc ] 


Vc 

J/4> 

Xo 


1.562(5) 
1.600(6) 
1.86(4) 


.645(1) 
.677(1) 
.686(3) 


.413(2) 
.424(2) 
.369(7) 


1.224(6) 
1.258(6) 
1.27(1) 


1.206(6) 
1.240(6) 
1.25(1) 
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fit to r 3 


m res 


m c 


fh c /m res 






Vc 
Xo 
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